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We present a coherent population trapping based scheme to attain sub-nanoscale resolution for atom local- 
ization, microscopy and lithography. Our method uses three-level atoms coupled to amplitude modulated probe 
field and spatially dependent drive field. The modulation of the probe field allows us to tap into the steep 
dispersion normally associated with electromagnetically induced transparency and offers an avenue to attain 
sub-nanometer resolution using just optical fields. We illustrate application of the techniques to the area of 
microscopy and lithography and show how multilevel schemes offer the possibility of improving resolution 
further. 



PACS numbers: 32.50.+d, 42.50.Gy, 32.80.-t 

Introduction. — The Rayleigh limit of resolution (A/2), 
when using a light field of wavelength A for imaging or for 
lithographic fabrication, primarily arises due to the choice of 
the optical design that is limited by diffraction. It has been 
shown that the ultimate limit of resolution is the Heisenberg 
limit and within its bounds, in principle, one can obtain sub- 
Ray leigh resolution [1]. In general, sub-Ray leigh resolution 
can be achieved by employing diffraction-free techniques. 
Several different avenues are being explored on widely dif- 
ferent fronts for obtaining sub-Rayleigh resolution based on 
near-field optics [2], quantum entanglement [3], monitoring 
fluorescence from atoms or molecules [4], and atomic co- 
herence effects [5-8]. Nanoscale spatial resolution is useful 
for applications to the fields of lithography [9], imaging [10], 
atom localization [11], high precision interferometry [12] and 
it offers precise spatial selectivity for atomic qubits necessary 
for quantum information processing [13]. 

We have shown recently that coherent population trap- 
ping (CPT) [7] can be used to localize atoms to sub- 
wavelength regimes [6]. The idea has been taken further 
by others to theoretically propose schemes for subwave- 
length microscopy [14], subwavelength patterning of Bose- 
Einstein condensates [15], nanoscale trapping potentials for 
atoms [16], and two dimensional localization by coupling the 
atom with two spatially dependent fields [17]. While the 
CPT based super-resolution techniques continue to be devel- 
oped; we note that a popular technique that attains nanome- 
ter resolution is the confocal fluorescence microscopy that 
makes use of the saturation of a two level transition. In an- 
other well-known, Stimulated Emission Depletion (STED), 
method one cuts down emission from outer regions by us- 
ing strong fields [4]. Fujita et al. [18] use the nonlinear re- 
lation between the excitation and fluorescence to demonstrate 
improved spatial resolution in three dimensions beyond the 
diffraction limit. 

Thus the nanoscale resolution is becoming rather common 
in microscopic and other applications. The question is how 
to break this barrier. In this letter we show how CPT can be 



adopted to reach sub-nanoscale resolution. We use the physics 
behind CPT and Electromagnetically-Induced Transparency 
(EIT) [8]. The sharp dispersion offered by EIT allows un- 
precedented control over the group velocity of the probe field. 
To take advantage of this steep dispersion we propose the use 
of amplitude modulated probe field and show that, in the best 
case scenario, the resolution can be further improved by a fac- 
tor of 20 beyond that offered by the unmodulated CPT. The 
modulated CPT gives resolution of 0.25 nm for the drive field 
of optical wavelength 500 nm. Thus we offer an avenue to 
attain sub-nanometer scale resolution which would have ap- 
plications to atom localization, imaging and lithography. 

The paper is organized as follows. First we briefly summa- 
rize the CPT-based atom localization scheme that is studied 
in detail in Ref. [6]. Then we introduce amplitude modula- 
tion for the weak probe field either as a perturbation or as a 
full modulation and study its effects on the localization. We 
present our results in terms of the point spread function of our 
scheme, which signifies the spread in the size of a point object 
as perceived by the scheme. 

Model. — We consider a three-level atom, as depicted in 
Fig. 1(a), illuminated by two optical fields, strong drive field 
on the 1 1)-|2) transition and weak probe field on 1 1)-|3) transi- 
tion with Rabi frequencies £l s and £l p respectively. We allow 
the fields to be detuned from their respective transitions but 
always assume two-photon resonance such that the two detun- 
ings are equal to a chosen value A. We also allow for decay of 
the excited state 1 1) into the other two states with decay rates 
y s and J p respectively. The Hamiltonian for this system is 

M 3 = -h(£l p |3) (1| + 0,(jc) |2) <l|)e- iA/ +H. c. . (1) 
The corresponding density matrix equations can be written as 

P = \W,p\- I ^(|l)(l|p-2p„|0|/)+p|l)(l|), 

n /=2,3 Z 

(2) 

with Yn = Ys and 713 = y p . In principle, the steady state solu- 
tions of the density matrix equations (2) can be obtained for 
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where 




FIG. 1 : (Color online) (a) The atomic levelscheme for coherent popu- 
lation trapping, (b) Comparison of the measurable feature size within 
the Rayleigh limit (white circle) and with the CPT scheme (small red 
circle in the center) for & = 1000 in the x-y plane. At this value of 
& the improvement in resolution is by a factor of about 50. 



the given Hamiltonian (1). Once the steady state is reached, 
such a system demonstrates coherent population trapping sig- 
nified by the form of the steady state given by = (£l p |2) — 
Cl s (x) |3))/Q, where £1 = y/\£l p \ 2 + \Q s (x)\ 2 ', this state does 
not evolve dynamically as Jtf =0. Thus, if the atom is 
initially prepared in the state |3), it will end up in the state 
at the steady state, as long as the two-photon resonance 
is maintained. As the population in the state can not es- 
cape, it is termed as the trapping state and the phenomena is 
called coherent population trapping. As shown in our earlier 
paper [6], taking the drive field to be a standing wave field, 
£l s (x) = £l s sinkx, the atomic dynamics becomes position de- 
pendent and possesses smaller-than-wavelength (A = 2n/k) 
features at the nodes of the standing wave field as witnessed 
by monitoring the population of the level |2): 



P22 



l + (£l s /£l p ) 2 l+^sin 2 **" 



(3) 



This function has the same structure as the transmission func- 
tion of the Fabry-Perot cavity [19], where the ratio, £& — 
|Q,s| 2 /|Qp| 2 , of the effective field intensities plays the role of 
the cavity finesse. The full width at half maximum (FWHM) 
of the curve is given by kAx = 2j\[M. In fig. 1 (b) we com- 
pare the feature-size measurable by the Rayleigh limit (A / 2) 
with the one measurable by the 2D version (see Fig. 5(a) 
for the corresponding levelscheme) of the CPT based scheme 
(X/nVM). For a moderate value of St = 1000, CPT-based 
scheme offers resolution improvement by a factor of 50. Fur- 
ther improvement is possible by increasing the value of 8%. 

Full Amplitude Modulation. — Now we consider full modu- 
lation of the probe field with simple standing- wave coupling 
field and once again monitor the population of level 1 2) , that is 
P22 as a measure of the point sperad function. The profiles of 
the coupling and modulated probe fields can be taken as £l s (x) 
and Clp(x,t) = £l p (x) cos vt. Thus, the population of the state 
1 2) can be determined to be 
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Some of the lower order terms are given by 
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i + [a p (x)/a s (x)] 2 



h = P22 



(i) _ 2 + 4[fl J (x)/fl / ,(x) j A ^ rA lrx 



l + [Q. p (x)/Q. s (x) 



-4[a s (x)/£l p (x)]\ 



- _ (2) _ 2[1 + 8 [Q 5 (x)/Cl p (x)] 2 (1 + [£l 5 (x)/£l p (x)] 2 ] 
U-P22 ~ / = 

+ 8[^(x)/^(x)] 2 (l + 2[^(x)/^(x)] 2 ). (6) 

The functions fu, which are functions of the transverse co- 
ordinates (x and also of y in the 2D case), are identified as 
the point spread functions of the proposed scheme for super- 
resolution. The FWHM of the point sperad function, which is 
a measure of how a point object is perceived by the method 
employed for microscopic detection, is a good measure of the 
resolution. We now consider the case where the drive field is a 
standing-wave, £l s (x) = £l s sin(fcjc), and the probe field has no 
position dependence Q. p (x) = £l p . The corresponding point 
spread functions (6) are plotted in Fig. 2 (a), with the general 
result that the higher-order Fourier components offer better 
resolution. The scheme can be extended to two dimensions 
(see Fig. 5 (a) for the corresponding level scheme) to obtain 
the plot of the FWHM of the corresponding point spread func- 
tions in 2D as in Fig. 2 (b). Taking the Fourier transform of the 




-0.06-0.04-0.02 0.00 0.02 0.04 0.06 

kx 



0.04-0.02 0.00 0.02 0.04 

kx 



FIG. 2: (Color online) Fully modulated CPT results with standing 
wave drive profile, (a) Point spread functions ft, Eq. (6), vs kx and 
(b) 2D FWHMs of the point spread functions against the background 
of the density plot of the 2D standing wave (sm 2 (kx) + sin 2 (ky)) 
near the node. Legend: Red (Solid): unmodulated CPT Result, Blue 
(Dashed): DC term (/ ), Green (Dotted): coefficient of cos2w (/ 2 ) 
and Black (Dot-Dashed): coefficient of cos4w (j^). The green and 
black plots are shifted vertically and may be inverted so that the 
FWHM can be visualized easily. Here & is larger by a factor of 
2 than the one used in the unmodulated CPT case due to time aver- 
aging required for the probe field, that is & = 2000. The resolution 
improvement from the Rayleigh limit of X/2 to the black circle is 
by a factor of 1000 and the resolution improvement from the red cir- 
cle (unmodulated CPT result to the black circle (coefficient of the 
second term in the Fourier expansion) is by a factor of about 23. 



(4) measured fluorescence signal from level 1 2) one can obtain the 
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FIG. 3: (Color online) Fully modulated CPT Results for transverse Laguerre-Gaussian (LG) drive profile, (a) With transverse drive and probe 
profile (b) zoom-in on (a), and (c) 2D version with the density plot of LG beam in the background. Legend: Red (Solid): unmodulated CPT 
Result, Blue (Dashed): DC term (/q), Green (Dotted): coefficient of cos2w (J2) an d Black (Dot-Dashed): coefficient of cos4w (fa), and 
Black (Thin Solid): guide to judge the FWHM. The blue and green plots are shifted vertically for easy comparison of FWHMs of different 
curves. Here we have chosen & = 1000. It is clear that the dot-dashed black curve corresponding to has the sharpest feature and offers 
the resolution improvement of about 20 compared to the unmodulated CPT result. 



term of a certain order. To illustrate, the point spread function 
corresponding to the cos(4w), fa, term gives a factor of 1000 
improvement over the Rayleigh limit, meaning this scheme 
can offer resolution of about 0.25 nm. 

For certain applications such as microscopy spatially mod- 
ulated drive and probe fields are more suitable in order to tar- 
get a smaller area of the sample to begin with. Therefore, 
in place of the standing wave modulation for the drive field 
we propose to use a donut (Laguerre-Gaussian) profile for the 
drive field and a gaussian profile for the probe beam. These 
profiles used on top of the perturbative modulation can lead 
to sharpening of the point spread function as discussed be- 
low. The spatial profiles can be written as: £l p (r) = e _r / 2w o, 

£l s (r) = \[M yjx 2 1 Wq e~ rl / 2w o , where wo is the beam waist 
and r is the transverse coordinate. The corresponding point 
spread functions (6) are plotted in Fig. 3 (a) and (b) along 



with the beam profiles in one dimension. It is clear that the 
spatial profiles offered by CPT are much finer than the spatial 
modulation of the light beams. This scheme can be easily ex- 
tended to two dimensions in a straightforward manner as the 
transverse beam profiles are naturally two dimensional. In the 
part (c) of the figure Fig. 3 we plot the FWHM as the size 
of the point spread function with the background showing the 
contour plot of the drive beam profile. The FWHM of the 
beam profile is too large to fit in the scale shown. When one 
monitors the second order term, in comparison with the un- 
modulated CPT results the resolution can be improved by a 
factor of about 3 and in comparison with the Rayleigh limit 
by a factor of about 150 for the chosen value of £% = 1000. 
Higher values of which is the ratio of the intensities of the 
drive and probe fields offers better resolution. 



Perturbative Amplitude Modulation. — We conclude the pa- 
per by examining a more traditional method of performing 
modulation spectroscopy. In this case the probe field is am- 
plitude modulated with a perturbative harmonic component 
sin vt where v is the frequency of the harmonic term. Whereas 
the drive field is spatially modulated along the x direction such 
that it takes the familiar form of a standing wave field to intro- 
duce position dependence in the atomic dynamics. The field 
Rabi frequencies can be written as: 

D. s (x) = Q. S sinkx, Cl p (t) = £l p (1 +asin vt) , (a«l). (7) 



Using a power series expansion in a for the density matrix 

elements pij = p-^ + ap-^ +a 2 p-^ H , the steady state 

solution of the density matrix equations of motion can be ob- 
tained. Once again the population of state |2) turns out to be 
position dependent and is given by: 
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(a 5 /n p ) 2 (3-(ci 5 /n p ) 2 
(i + (iyo p )2) 3 



(asinvt) 2 



(8) 



In Fig. 4 we plot the coefficients enclosed in the square brack- 
ets in the above relation, for the field profiles given in Eq. (7), 
versus the dimensionless position coordinate kx. For the cur- 
rent system the point spread function is proportional to the 
appropriate coefficients determined above. It is clear that the 
resolution is improved ever so slightly by employing a higher 
and higher order term in the expansion of population. In an 
experimental setting, P22 will be measured as a whole. Then 
knowing the particular values used for the perturbation param- 
eter a and the harmonic frequency v with the help of curve- 
fitting the required term can be determined by substituting for 
the lower order terms with their theoretical values. The per- 
turbative amplitude modulation can also be applied to the case 
where drive field has the Laguerre-Gaussian beam profile and 
the probe is Gaussian. Similar results to the case of full mod- 
ulation are obtained that offer finer resolution than the trans- 
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we take advantage of the dispersion accompanying sharp res- 
onances in EIT that is closely related to the phenomena of 
CPT. The improvement in resolution obtained in comparison 
with the Rayligh limit is by a factor of about 1000 in the best 
case scenario discussed here which gives sub-nanometer res- 
olution with the use of optical fields. 



FIG. 4: (Color online) Perturbatively modulated CPT Results, (a) 
ID, and (b) 2D version with the density plot of the 2D standing wave 
near the node. Legend: Red (solid): Unmodulated CPT Result or 

p$, Blue (Dashed): p^ } , Green (Dotted): p^\\ and Black (Thin 
Solid): guide to judge the FWHM. The blue (dashed) and green (dot- 
ted) plots are shifted vertically for easy comparison and visualization 
of FWHM. Here we have chosen & = 1000. The resolution improve- 
ment for the first order (Green) is by a factor of 3 compared with 
unmodulated CPT result in Red. 



verse spatial profile of the beams. 

New possibilities. — The essentially ID model presented in 
Fig. 1 can be extended to 2D by taking advantage of the mul- 
tilevel nature of atoms as shown in Fig. 5(a). Two standing 
waves in x and y directions induce sharp features in P33 al- 
lowing sub-nanoscale localization in the x-y plane. Further, 
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FIG. 5: (a) Multilevel scheme for 2D imaging. The population of 
level 3 is P33(x,y) and shows sharp features. Notice the x and y 
dependance of the drive field coupling to level |3) (b) Multilevel 
scheme for further enhancement of resolution in a single dimension. 
This can be envisioned to be two coupled A systems giving square of 
the Fabry-Perot function for p\\(x). 

considering the lower level states |2) and |3) as hyperfine 
sublevels of a real atom with F > 2 one can obtain multiple 
A systems in a single atom that are simultaneously excited; 
such systems have been shown to be useful for generation 
of skyrmions in the context of BEC and Laguerre-Gaussian 
beam interaction [20]. In the case shown in Fig. 5(b) the re- 
sult with unmodulated CPT gets modified to a function that 
can be approximated by the square of the Fabry-Perot trans- 
mission function (3) and hence is sharper by a factor of about 
1.5 in the case of unmodulated CPT. Adding the probe field 
modulation on top of this CPT in the coupled A system one 
can, in principle, improve the resolution further. The ideas 
presented so far can also be extended to three dimensions by 
employing standing waves in the three directions much on the 
same lines as techniques for optical trapping of neutral atoms 
in 3D [21]. 

Conclusions. — We offer sub-nanoscale resolution by mod- 
ulating probe fields in standard CPT. By modulating the CPT 
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